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ABSTRACT 



Recently, we have developed a method useful for mapping large-scale horizontal velocity fields in the solar photosphere. The method 
was developed, tuned and calibrated using the synthetic data. Now, we applied the method to the series of Michelson Doppler Imager 
(MDI) dopplergrams covering almost one solar cycle in order to get the information about the long-term behaviour of surface flows. 
We have found that our method clearly reproduces the widely accepted properties of mean flow field components, such as torsional 
oscillations and a pattern of meridional circulation. We also performed a periodic analysis, however due to the data series length and 
large gaps we did not detect any significant periods. The relation between the magnetic activity influencing the mean zonal motion 
is studied. We found an evidence that the emergence of compact magnetic regions locally accelerates the rotation of supergranular 
pattern in their vicinity and that the presence of magnetic fields generally decelerates the rotation in the equatorial region. Our results 
show that active regions in the equatorial region emerge exhibiting a constant velocity (faster by 60 ± 9 ms~' than Carrington rate) 
suggesting that they emerge from the base of the surface radial shear at 0.95 R , disconnect from their magnetic roots, and slow down 
during their evolution. 

Key words. Sun: photosphere - Sun: magnetic fields - Sun: activity 



1. Introduction 

The largest scale velocity fields in the solar photosphere consist 
of rotation profile and a meridional flow pattern. 

Basically, the differential rotation is described as an integral 
of the zonal component of the studied flow field. The inte- 
grated flow field may be obtained using spectroscopic method, 
using tracer-type measurements, or using helioseismic inver- 
sions. The latest case allows to measure the solar rotation not 
only as a function of heliographic latitude, but also as a func- 
tion of depth. From the helioseismic inversion we know that 
throughout the convective envelope, the rotation rate decreases 
monotonically toward the poles by about 30 %. Angular veloc- 
ity contours at mid-latitudes are nearly radial. Near the surface 
at the top of the convection zone there is a layer of a large radial 
shear in the angular velocity. At low and mid-latitudes there is an 
increase in the rotation rate immediately below the photosphere 
which persists down to r ~ 0.95 R Q . The angular velocity varia- 
tion across this layer is roughly 3 % of the mean rotation rate and 
accord ing to the helioseismic analysis of ICorbard & Thompsonl 
(I2002I) the angular velocity a> decreases within this layer approx- 
imately as r , where r is a radial coordinate. At higher latitudes, 
the situation is less clear. For the overview of solar differential 
rotation measurements see ISchroterl(ll985l) or a more recent re- 
view bv lBeckl(l2000l) . 

The torsional oscillations, in which narrow bands of 
faster than average rotation, interpreted as zonal flows, mi- 
grate towards the solar equator during the sunspot cycle, 
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were discovered bv | Howard & Labontd (Il980h . Later research 
(lAntia & Basul |200lT found that there exist two different 
branches of torsional oscillations. At latitudes below about 40 °, 
the bands propagate equatorward, but at higher latitudes they 
propagate poleward. The low-latitude bands are about 15 wide 
in latitude. The flows w ere studied in surface Doppler mea- 
surements dUlrichl 120011). an d also using local helioseismology 
dKosovichev & Schoul 19971) . The surface pattern of torsional 
oscillations penetrate deep ly in the convec t ion zo ne, possibly 
to its base, as suggested bv IVorontsov et all (120021) . The ampli- 
tude of the angular velocity variation is about 2-5 nHz, which is 
roughly 1 % of the mean rotational rate (5-10 m s _1 ). The direct 
comparison betwee n different t e chniq ues inferring the surface 
zonal flow pattern dHowe et al 1 12006b showed that the results 
are sufficiently coherent. 

The surface magnetic activity corresponds well with the tor- 
sional oscillation pattern - the magnetic activity belt tends to lie 
on the poleward side of the faster-rotating low-latitude bands. 
The magnetic activity migrates towards the equator with the low- 
latitude bands of the torsional osci l lation s as the sunspot cy- 
cle progres ses (IZhao & Kosovichevl 12001 . Some studies (e. g. 
iBeck et aTll2002h suggest that meridional flows may diverge out 
from the activity belts, with the equatorward and poleward flows 
well correlated with the faster and slower bands of torsio nal os- 
cillations. In a recent theoretical study bv lRempel d2007l) it was 
suggested that the poleward-propagating high-latitude branch 
of the torsional oscillations can be explained as a response of 
the coupled differential rotation/meridional flow system to peri- 
odic forcing in midlatitudes of either mechanical (Lorentz force) 
or thermal nature. The equatorward-propagating low-latitude 
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branch is most likely not a consequence of the mechanical forc- 
ing alone, but rather of thermal origin. 

The axisymmetric flow in the meridional plane is generally 
known as the meridional circulation. The meridional circulation 
in the solar envelope is much weaker than the differential rota- 
tion, making it relatively difficult to measure. 

Two principal methods are widely used to measure the 
meridional flow: feature tracking and direct Doppler measure- 
ment. There are several difficulties complicating the measure- 
ments of the meridional flow using tracers. Sunspots and fila- 
ments do not provide sufficient temporal and spatial resolution 
for such studies. Sunspots also cover just low latitudinal belts 
and do not provide any informations about the flow in higher 
latitudes. Doppler measurements do not suffer from the problem 
associated with the tracer-type measurements, however they in- 
troduce another type of noisy phenomena. It is difficult to sepa- 
rate the meridional flow signal from the variation of the Doppler 
velocity from the disc centre to the limb. Using different tech- 
niques, the parameters of the meridional flow show large dis- 
crepancies. It is generally assumed that the solar meridional flow 
in the close subphotospherical layers is directed poleward with 
one cell per hemisphere. Such flow is a lso produced by early 
global hydrodynamical s imulation such as l Glatzmaier & Gilmanl 
dl982l) . As reviewed by lHafhawavl d 19961) . the surface or near 
sub-surface velocities of the meridional flow are generally in the 
range 1-100 ms _1 , the most often measured values lie within 
the range of 10-20 ms~'. The flow has often a complex lat- 
itudinal structure with both poleward and equatorward flows, 
mu ltiple cells, and large asym metries with respect to the equa- 
tor. IZhao^^^soviche^ (|2003) used the time-distance helioseis- 
mology to infer the properties of the meridional flow in years 
1996-2002. They found the meridional flows of an order of 
20 ms~', which remained poleward during the whole period of 
observations. In addition to the poleward meridional flows ob- 
served at the solar minimum, extra meridional circulation cells 
of flows converging toward the activity belts are found in both 
hemispheres, which may imply plasma downdrafts in the activ- 
ity belts. These converging flow cells migrate toward the so- 
lar equa tor together with the activity belts as the solar cycle 
evolves. iBeck et aT] d2002l) measured the meridional flow (and 
torsional oscillations) using the time-distance helioseismology 
and found the residual meridional flow showing divergent flow 
patterns around the solar activity belts below a depth of 18 Mm. 

The most complete maps of the torsional oscillations and the 
meridional flow available at present have been constructed on the 
basis of Mt. Wilson daily magnetograms (see lUlrich & Bovdenl 
120051) . The measurements cover more than 20 years (since 1986) 
and the results obtained using this very homogenous material 
agree well with the properties described above. 

The modern dynamo flux-transport models use the merid- 
ional flow and the differential rotation as the observational in- 
put. In the models by Dikpati et al. dDikpati et al. 1 120061 or 
iDikpati & Gilm an 2006) the return meridional flow at the base 
of the convection zone is calculated from the continuity equa- 
tion. They found the turnover time of the single meridional cell 
of 17-21 years. The meridional flow is assumed to be essential 
for the dynamo action, global magnetic field reversal and fore- 
cast of the future solar cycles. 

There are known many relations of the differential rotation 
profile to the ph as e of the progre s sing s olar cycle - see e.g. 
iJavaraiaril d2003l) or lJavaraiah et alj d2005l) - showing for exam- 
ple different properties of the differential rotation profile in the 
odd and even solar cycles. The rotation of the sun spots in relation 
to their morphological type was studied e. g. by iBalthasar et alj 



d!986l) who found that more evolved types of su nspots (E, F, G 
and H type) rotate slower than less evolved types. iRuzdiak et aT] 
d2004l) investigated Greenwich Photoheliographic Results for 
the years 1874-1976 and found a clear evidence for the decel- 
eration of th e sunsp ots in the photosphere with their evolution. 
iHerdiwijaval d2002l) found that the leading part of a complex 
sunspot group rotate about 3 % faster than the following part. 
The dependence of the rotation of sunspot on their size and posi- 
tion in the bipolar region was investigated bv lD'Silva & Howard! 
d 19941) . They explained the observed behaviour through a sub- 
tle interplay between the forces of magnetic buoyancy and drag, 
coupled with the role of the Coriolis force acting on rising flux 
tubes. This dynamics of rising flux tubes also explains the faster 
rotation of smaller sunspots. In average, sunspots rotate about 
5 % faster than the surroun ding plasma . 

In the theoretical study (Brun, 20041) based on 3-D numerical 
simulations of compressible convection under the influence of 
rotation and magnetic fields in spherical shells, the author stated 
that in the presence of magnetic field the Maxwell stresses may 
oppose the Reynolds stresses and therefore the angular momen- 
tum is propagated more to the poles than without the presence 
of magnetic fields. As a consequence the rotation profile is more 
differential in the periods of lowered magnetic activity and it 
leads to the increase of the rotation in lo w latitudes. This be- 
haviou r was observed in many studies, e. g. Hathaway & Wilsonl 
<f!990h . 

The subject of this work is a verification of the perfor- 
mance of the method described in lSvanda et al.l d2006l) (hereafter 
Paper I) on the real data and the investigation of long-term prop- 
erties of the flows at largest scales obtained with this method. 
We shall also discuss the influence of magnetic fields on the 
measured zonal flow in the equatorial region. This topic will be 
studied more in detail in one of the next papers in the series. 

2. Data processing 

The horizontal photospheric velocity fields a re calculated using 
local correlation tracking (LCT) algorithm dNovemberl |1986[) 
applied to series of processed full-dis c dopplergrams meas ured 
by Michelson Doppler Imager (MDI; IScherrer et all 1 1995b on- 
board Solar and Heliospheric Observatory (SOHO). In the dopp- 
lergrams, the supergranular pattern is tracked in order to obtain 
properties of the velocity field of the larger scale. 

For this study we have processed all suitable data measured 
by MDI. The instrument observed full-disc dopplergrams ap- 
proximately two months each year in a high cadence of one 
frame per minute. These Dynamic campaigns provide suitable 
material for our method. Between May 23, 1996 and May 22, 
2006 we have 806 days covered by high-cadence measurements. 
In each of these days, two 24-hour averages sampled every 12 
hours were calculated. In some days, MDI had significant gaps 
in measurements, so in such cases we did not have enough ho- 
mogeneous material to process. Therefore, in all the Dynamic 
campaigns 502 days were useful for our analysis and we cal- 
culated 1004 full-disc horizontal velocity fields. The processing 
of almost 3 TB of primary data took several months using fast 
computers running in the network of W. W. Hansen Laboratory, 
Stanford university. 

The data were processed using the technique described in de- 
tail in Paper I. We bring just a brief summary. The method pro- 
cesses 24-hours series of MDI full-disc dopplergrams, contain- 
ing 1 440 frames. The one-day series first undergoes the noise 
and disturbing effects removal. From all frames, the line-of-sight 
component of the Carrington rotation is subtracted and the effect 
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Fig. 1. Mean meridional flow in time and heliographic latitude. It can be clearly seen that for almost all the processed measurements 
a simple model of one meridional cell per hemisphere would be sufficient. However, some local corruptions of this simple idea can 
be noticed on both hemispheres. 
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Fig. 2. Torsional oscillations. The residua of the mean zonal flow with respect to its parabolic fit displayed in time and heliographic 
latitude. In the period of weak magnetic activity the pattern of belts propagating towards the equator is very clear. In the periods of 
stronger magnetic activity the flow field is influenced by the local motions in active regions and therefore the pattern of torsional 
oscillations is not clearly seen. 



of a perspective is corrected. The frames are transformed so that 
the heliographic latitude of the disc centre bo = and the posi- 
tion angle of the solar rotation axis P = 0. Then, the p-modes of 
the solar oscillati ons are removed using a weighted average (see 
lHathawavLfl988h . The weights have a Gaussian form given by 
the formula: 

(Art 2 I b 2 - (At) 2 \ 

w(At)=e^ ^a2 )' W 

where Af is a time distance of a given frame from the central one 
(in minutes), b — 16 minutes and a — 8 minutes. We sample 
averaged images in the interval of 15 minutes. The filter sup- 
presses more than five hundred times the solar oscillations in the 
2-4 mHz frequency band. 



The processing of averaged frames consists of two main 
steps. In the first main step the mean zonal velocities are cal- 
culated and, on the basis of expansion to the Fay's formula 
lj — co + c\ sin 2 b + C2 sin 4 b, the differential rotation is removed. 
In the second main step, the LCT algorithm with an enhanced 
sensitivity is applied. Finally, the differential rotation (obtained 
in the first step) is added to the vector velocity field obtained in 
the second main step. Both main steps can be divided into sev- 
eral sub-steps, which are mostly common. 

1. The data series containing 96 averaged frames is "derotated" 
using the Carrington rotation rate in the first step and using 
the calculated differential rotation in the second step. 

2. Derotated data are transformed into the Sanson-Flamsteed 
coordinate system to remove the geometrical distortion 
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caused by the projection to the disc. The Sanson-Flamsteed 
(also known as sinusoidal) pseudo-cylindrical projection 
conserves the areas and therefore is suitable for the prepa- 
ration of the data used by LCT. 

3. Rema pped data undergoes the k-co filtering (e. g. iTitle et al.L 
1 1 989b with the cut-off velocity 1 500 ms" 1 for suppression 
of the noise coming from the evolutionary changes of super- 
granules, of the numerical noise, and for the partial removal 
of the "blind spot" (an effect at the centre of the disc, where 
the supergranular structures are almost invisible in doppler- 
grams due to the prevailing horizontality of their internal ve- 
locity field). 

4. Finally, the LCT is applied: the lag between correlated 
frames is 4 hours, the correlation window has FWHM 60", 
the measure of correlation is the sum of absolute differences 
and the nine-point method for calculation of the subpixel 
value of displacement is used. The calculated velocity field 
is averaged over the period of one day. 

5. The resulting velocity field is corrected using the formula 

Vcoi = 1.13 Vcaic, (2) 

where v ca i c is the magnitude of velocities coming from the 
LCT, and v cor the corrected magnitude. The directions of the 
vectors before and after the correction are the same. The cal- 
ibration formula was obtained from the tests on the synthetic 
data (see Paper I). Finally, v x component is corrected for the 
data-processing bias of - 1 5 m s~ 1 determined in Paper I too. 

In this study, we were interested only in the properties of the 
mean zonal and meridional components. Therefore, from each 
two-component horizontal velocity field the mean zonal and 
meridional component depending only on heliographic latitude 
were calculated as the longitudinal average of the flow map, us- 
ing 135 longitudinal degrees around the central meridian. As 
stated in Paper I, the accuracy for each velocity vector is 15 m s~' 
for velocities under 100 m s and 25 m s~' for velocities above 
100 ms" 1 . These inaccuracies have a character of a random er- 
ror, therefore for the mean zonal and meridional components the 
accuracy is in the worst case 1 m s~' . 

Th e performance of the method was verified bv lSvanda et alj 
d2007l) . The results of comparisons between the technique used 
in the present study and the time-distance helioseismology show 
that both methods reasonably match. The calculated surface 
flows may be biased by the projection effects, although the tests 
on the synthetic data ( Paper I) did not show any signs of them. 
lHathawav et al.l d2006l) showed that the apparent superrotation 
of the structures tracked in dopplergrams reported by many stud- 
ies can be explained as the projection effect. However, this bias 
would produce the systematic error, which should influence nei- 
ther the periodic analysis, nor the relative motions of the active 
regions with respect to its surroundings. 

3. Results 

3. 1 . Long-term properties 

For the study of long-term evolution of surface flows maps con- 
taining the mean zonal and meridional components were cal- 
culated. The maps of mean meridional component in time and 
heliographic latitude are shown in Fig. Q] We can clearly see 
that on the northern hemisphere dominates the flow towards the 
northern pole while on the southern hemisphere the flow towards 
southern pole prevails. The "zero line", the boundary between 
the flow polarity is not located exactly on the solar equator and 



seems to be shifted to the south in the period o f increased so- 
lar act ivity (2001 and 2002). In agree ment with iMeunier et al.l 
(119971) andlC ameron & Hopkins! (1 1998) we found the meridional 
flow stronger in the periods of increased solar activity by about 
10ms" 1 than in periods with lower magnetic activity. 

A similar map was made in the same way for the zonal com- 
ponent. The mean equatorial zonal velocity for all the data is 
1900 ms~'. For all the processed data the dependence on lati- 
tude is close to a parabolic shape, parameters of which change 
slowly in time. The residua of the zonal velocity with respect to 
its parabolic fit given by 

Vb = ao + a\b + 02b 2 , (3) 

where b is the heliographic latitude and v\, the mean zonal veloc- 
ity in the given latitude, were calculated in order to see if we are 
able detect torsional oscillations in our measurements. As it is 
displayed in Fig. [2] the method clearly reveals torsional oscilla- 
tions as an excess of the mean zonal velocity with respect to the 
zonal velocity in the neighbourhood. The behaviour of torsional 
oscillations is in agreement with their usual description - the ex- 
cess in magnitude is in order of 10 m s _1 , they start at the begin- 
ning of the solar cycle in high latitudes and propagate towards 
the equator with the progress of the 1 1 -year cycle. However, with 
our method the visibility of torsional oscillations decrease with 
increasing solar activity. In the periods of strong activity both 
belts are not so clearly visible since the large-scale velocity field 
and its parabolic fit are strongly influenced by the presence of 
magnetic regions. However, the torsional oscillations belts still 
remain visible when the mean zonal component is symmetrised 
with respect to the solar equator. We did not focus on study of 
meridional flow or torsional oscillations depending on time and 
latitude, we just used them to check the ability and performance 
of our method. 

3.2. Periods in the mean components 

The mean zonal and meridional components in the equatorial 
area (averaged in the belt b = -5 - +5 °) were analysed in or- 
der to examine the periods contained in the data. Since the data 
are not equidistant at all, we cannot us e a simple harmonic anal- 
ysis, so that the Stellingwerf method ( Stellingwerf , 1978) was 
applied. It works on the principle of the phase dispersion min- 
imization. The method sorts the data for every searched period 
into the phase diagram. Then the phase diagram is divided in a 
few (mostly ten) parts and for every part the mean dispersion 
is calculated. If the studied period is reasonable, the data-points 
group along the periodic curve and the dispersion in each part 
of the phase diagram is smaller than the dispersion of the whole 
data series. The normalized parameter 9 e (0, 1) describes the 
quality of a given period. 

We cannot completely exclude the influence of the calculated 
flow fields by the the position and orientation of the solar disc 
(position angle of the rotation axis P, heliographic latitude of the 
centre bo), so we put also the series of both parameters sampled 
on the same dates when our measurements of surface flows exist 
through the period analysis. We also verified the periods caused 
by sampling of the data using the same method. Periodograms 
are displayed in Fig. [3] 

We conclude that we did not detect any significant period in 
the available data set. As it can be seen in Fig. [3] there exist non- 
convincing (the values of the parameter 6 are quite high, which 
means that the period probably is not significant) signs of periods 
detected in the real data, which are not present in the control data 



M. Svanda et al.: Large-scale horizontal flows in the solar photosphere II: Long-term behaviour. 



5 




0.6 



0.4 



0.2 



500 



1000 



1500 2000 
Period [days] 



2500 3000 



3500 




1500 2000 
Period [days] 



3500 




0.6 



0.4 



0.2 



0.6 



0.4 



0.2 



500 



1000 



1500 2000 
Period [days] 



2500 3000 



3500 



500 



1000 



1500 2000 
Period [days] 



2500 3000 



3500 



Fig. 3. Periodograms determined using the Stellingwerf method. Parameter 6 signifies the normalized phase variation. Upper left: 
Periodogram of mean equatorial zonal velocity. Upper right: Periodogram of mean equatorial meridional velocity. Bottom left: 
Periodogram of sampled heliographic latitude of the centre of the solar disc. Bottom right: Periodogram related to the sampling of 
data. 



set. Values of the suspicious periods are 657 days (1.80 years) 
in the meridional component and 1712 days (4.69 years) in the 
zonal co mponent. We note tha t the 1.8-year period was also de- 
tected by iKnaack et alj d2004f) . It is claimed to be related to a 
possible Rossby wave r-mode signature in the photosphere with 
azimuthal order m ~ 50 reported by iKuhn et ail d2000t) . but 
lately disputed e. g. by IWilliams et alj d2006l) . The period esti- 
mate for such an r -mode is close to 1.8 years. According to 
IKnaack et alj d2005l) . such a periodicity was observed in the total 
magnetic flux only on southern hemisphere from 1997 to 2003. 
The coupling between the zonal flow and the meridional circula- 
tion could transfer the signal of the r-mode motion to the mean 
meridional component. 

The detected suspicious periods may not be of solar origin. 
The sparse data set suffers from aliasing caused by a bad cov- 
erage of the studied interval. To confirm the periods far more 
homogenous data set is needed. This may be a task for ongoing 
space borne experiment Helioseismic Michelson Imager (HMI), 
which will be a succesor of MDI. 

Detected periodicities are absent in Mt. Wilson torsional os- 
cillations time series, which is far more homogenous ma terial 
than the one used in this study. ISnodgrass & Ulrichl d 19901) also 
did not find any time variations in the study tracking the features 
in the low resolution dopplergrams covering homogenously 20 
years of Mt. Wilson observations. All the arguments written 



above led us to leave the detected periods as suspicious, as they 
cannot be confirmed from the current data set. 



3.3. Relation to the magnetic activity 

We also investigated the coupling of equatorial zonal velocity 
(average equatorial solar rotation) and the solar activity in the 
near-equatorial area (belt of heliographic latitudes from -10° 
to +10°). The average equatorial zonal velocity incorporates 
the average supergranular network rotation and also the move- 
ment of degenerated supergranules influenced by a local mag- 
netic field with respect to their non-magnetic vicinity. Indexes of 
the solar activity were extracted from the daily reports made by 
Space Environment Center National Oceanic and Atmospheric 
Administration (SEC NOAA). Only the days when the measure- 
ments of horizontal flows exist were taken into account. As 
the index of the activity we have considered the total area of 
sunspots in the near-equatorial belt and also their type. 

First of all we computed the correlation coefficient p between 
the mean equatorial zonal velocity and the sunspot area in the 
near-equatorial belt and got a value of p = -0.17. We cannot 
conclude that there is any linear relation between these two in- 
dices. The dependence of both quantities is plotted in Fig. [4] 
We can clearly find two different regimes which are divided by 
the velocity of approximately 1890 ms _1 . In one regime (77 % 
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Fig. 5. Sunspot area sampled in the same times when the measurements of the horizontal flows exist. Two regimes of the near- 
equatorial belt rotation are displayed. Diamonds denote "fast" rotating equatorial belts, crosses the "scattered" group. 
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Fig. 4. Mean zonal equatorial velocity versus the sunspot area 
in near-equatorial belt. We decided to divide the data in two 
regimes along the velocity axis. Although the division is arbi- 
trary, we believe that it is supported by the theory of the dynam- 
ical disconnection of sunspots from their roots. 



of the cases), the equatorial belt rotates about 60 ms~' faster 
(1910 + 9ms~ 1 ; hereafter "fast group") than Carrington rotation, 
in the other one (23 %) the rotation rate is scattered around the 
Carrington rate (1860 + 20 ms~'; hereafter "scattered group"). 
The division in two suggested groups using the speed criterion 
is arbitrary. If there exist only two groups, they certainly overlap 
and only a very detail study could resolve their members. We 
may also see more than two groups in Fig. [4] The arguments for 
division in just two groups will follow. 

For both regimes, there does not exist a typical sunspot area. 
The distribution of both regimes in time is displayed in Fig. [5] 
The data in the periods of larger solar activity (years 2001 and 
2002, these are also the only years when the data cover two 
Carrington rotations continuously) show that both regimes al- 
ternate with a period of one Carrington rotation. 

The histogram of the mean zonal equatorial velocity has a 
similar, i. e. bimodal, character like in Fig. [4] with a greater sec- 
ond peak, because such a histogram is constructed not only from 



belts cointaining magnetic activity but also from the belts where 
not magnetic activity was detected. The mean equatorial rotation 
for all the data is 1900ms -1 , 1896ms" 1 for the equatorial rota- 
tion in a presence of sunspots, and 1904 ms" 1 for days without 
sunspots in the equatorial region. 

Such bimodal velocity distribution is in disagreement 
with the results o btained by time-distance helioseismology by 
IZhao et al] d2004f) . In this work authors found that the stronger 
the magnetic field the faster such a magnetic element rotates. 
They observed about 70 ms" 1 faster rotation than the average 
for magnetic areas with magnetic fields stronger than 600 G. It 
might be possible that the size of the magnetic area and its mag- 
netic field strength play di fferent ro l es in t he influencing of the 
plasma motions. However. iMeunie ri (120051) showed that the size 
of the magnetic area and the maximum magnetic field strength 
or the total magnetic flux correlate quite well. 

Detailed studies of the sunspot drawings obtained from the 
Patrol Service of Ondfejov Observatory and the Mt. Wilson 
Observatory drawings archive revealed that in the "fast" group, 
the new or growing young active regions were present in the 
equatorial belt. On the contrary, in the "scattered" group the 
decaying or recurrent active regions prevailed in the equato- 
rial area. The deceleratio n of the sunspot grou p with its evolu- 
tion was noticed e. g. by iRuzdiak et alj d2004l) . Moreover, our 
results suggest that the new and rapidly growing sunspots in 
the studied sample (March to May 2001 and April to June 
2002) move with the same velocity. This behaviour could be ex- 
plained by emergence of the local magnetic field from a con- 
fined subphotospheric layer. According to the rough estimate 
dCorbard & Thompsoni [20021) the speed of 1910 + 9 ms" 1 corre- 
sponds to the layer of 0.946 ±0.008 R Q , where the angular veloc- 
ity of rotation suddenly changes. During the evolution, the mag- 
netic field is disrupted by the convective motions. An interesting 
behaviour is displayed by the alternation of the "fast" and "scat- 
tered" regimes (see Fig. [6} with the period of one Carrington 
rotation. It suggests that active regions in the equatorial region 
emerge in groups. We have to keep in mind that our study pro- 
duces a very rough information due to the averaging of all efects 
in the equatorial belt. 

The observed behaviour could be a manifestation of the 
disconnection of magnetic field lines from the base of the 
surface shear during evolution of the gro wing sunspot group. 
This b ehaviour was theoretically studied by lSchussler & Rempell 
d2005l) . They suggested the dynamical disconnection of bipolar 
sunspot groups from their magnetic roots deep in the convection 
zone by upflow motions within three days after the emergence 
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Fig. 6. Left: Distribution of two equatorial rotation modes in year 2002. Right: Derivatives of the mean zonal velocity (solid curve) 
and the sunspot area in the near-equatorial region (dashed curve) in year 2002. Both quantities correlate with each other quite nicely. 
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Fig. 7. Active region morphological types (left) and the number of active regions in the near-equatorial belt distributed as a function 
of time and sunspot area. 



of the new sunspot group. The motion of sunspots changes dur- 
ing those three days from "active" to "pasive". The active mode 
is displayed by motions reasonable faster with respect to a non- 
magnetic origin. The passive mode means mostly the deceler- 
ation of sunspot motions and influence of the sunspot motions 
only by the shallow surface plasma dynamics. The theory of the 
disconnection of sunspot groups from their magnetic root sup- 
ports the division of the data set in two groups. 

As an example, we selected the active region NOAA 9368 
(Fig. [8]l to show the behaviour of large-scale velocities in time. 
We see that the leading part of the active region rotates faster 
than the surroundings in the first day of observation and the 
whole group slows down in next two days. The inspection of 
details in the behaviour of selected active regions in the whole 
data set will be the subject of the ongoing studies. 

Division of the equatorial belt into 10 sectors helped the 
investigation of the behaviour of different active region types. 
The mean zonal velocity in the sector containing the studied 
active region according to morphological type of active re- 
gions are summarized in Table \T\ It can be clearly seen that, 
in average, more evolved active regions rotate slower than less 
evolved or young acti ve regions, what is in agreement with e. g. 
Balt hasar et al.l d 1 986b . We have also replotted Fig. |4]using seg- 
mented equatorial belt and different active region types and got a 
very similar result. We did not find any particular behaviour for 
various active region types (see Fig. [7j>. 

We have also focused how the presence of the magnetic ac- 
tive areas will influence the average flow field. Since we found 
that a direct correlation is weak due to existence of two different 
regimes, we decided to study the temporal change of both quan- 
tities. The aim is to study whether an emerging active region in 
the near-equatorial belt will influence the average equatorial ro- 
tation. We computed numerical derivatives of the total sunspot 
area in the near-equatorial belt and of the average zonal equato- 
rial flow. We have found that the correlation coefficient between 
both data series is p — 0.36 and is higher for the "fast group" 



(p — 0.41) than for the "scattered group" (p = 0.24). The corre- 
lation is higher in periods of increased magnetic activity in the 
equatorial belt. For example, for data in year 2001 the correla- 
tion coefficient is P2001 = 0.58 and for year 2002 P2002 = 0.52; 
see Fig. [6] In both particular cases, the correlation is higher for 
the "fast regime" (p ~ 0.7) than for the second group. 

It is important to know that LCT in our method measures ba- 
sically the motions of supergranules influenced by the magnetic 
field, not by spots that are recorded in the used solar activity in- 
dex. Therefore the results can be biased by the fact that the pres- 
ence of magnetic field does not necessarily mean the presence of 
sunspot in the photosphere. 

We think that despite an apparent disagreement, our ap- 
parently conflicting results can be valid and are in agreement 
with the results published e arlier. Basically, as described b y e. g. 
iHathawav & Wilson! ( I 19901) and explained in the model of lBrunl 
(12004 . the solar rotation in the lower latitudes is slower in 
the presence of magnetic field. This is summarized in Table Q] 
and displayed in Fig. |4] In most cases, "spotty" equatorial belts 



Table 1. Mean synodic rotation velocities of different active re- 
gions types and the average equatorial rotation of all data. The 
measured speed v alues may be systemati cally biased by the pro- 
jection effect (see lHathawav et alll2006l) . 



Sunspot 


Mean rotation 


type 


[m s" 1 ] 


A 


1893 ±48 


B 


1893 ±49 


C 


1890 ±71 


D 


1880 ±73 


E 


1880 ±50 


F 


1874 ±73 


H 


1872 ±51 


Average 


1900 
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Fig. 8. Case study: Evolution of the flows in and around the ac- 
tive region NOAA 9368 on March 6, March 7, and March 8, 
2001 . The leading polarity rotate significantly faster than the fol- 
lowing one and the non-magnetic surroundings in the first day. 
The whole group slows down in the other two days. As the back- 
ground image, the MDI magnetogram smoothed to the resolu- 
tion of the measured flow field is used. The magnetic field in- 
tensities are displayed in the range of -1800 (black) to +1800 
(white) Gauss in the linear scale. 



seem to rotate slower than the average for the whole data series. 
However, it is clear that emerging active regions cause in most 
cases the increase of the rotation rate. This is in agreement with 
a generally accepted statement found first bv lHoward & Harvevl 



d 19701) and lGolub & Vaianaldl978l) . The relation, obtained using 
linear fit on our data set, can be described by the equation 

Av ~ 0.2 AA sunspots m s" 1 , (4) 

where Av is a change of the equatorial rotation speed with re- 
spect to the Carrington rotation and AA sunspots is a relative change 
of sunspot area in the equatorial belt (in 10~ 6 of solar hemi- 
sphere). We estimate that strong local magnetic areas rotate few 
tens of ms~' faster than the non-magnetic surroundings. 

The difference in the behaviour of the flow fields in the re- 
gions occupied by the magnetic fi eld and their vicinity was stud- 
ied also e. g. bv lTitle et alj ( 1 19921) using the high-resolution data 
obtained at the Swedish Vacuum Solar Telescope on La Palma, 
Canary Islands. The authors found that the magnitude of hori- 
zontal velocities measured by LCT on the granular scale is in 
the regions of the quiet Sun larger than in an active-region plage. 
The high-resolution velocity fields are of different nature than 
the large-scale ones studied by our method. Flow fields on gran- 
ular scales are mostly chaotic due to the turbulent behaviour. In 
the regions occupied by the magnetic field, the motions become 
more organized, the chaotic component is suppressed, and there- 
for e the amplitude of the horiz ontal velocity is generally lower. 

ISnodgrass & Ulrichl (1 19901) tracked the features in the low- 
resolution dopplergrams measured at the Mt. Wilson 46 m tower 
telescope for the period of 20 years. They interpreted the de- 
tected flows as the velocity of the supergranular network, al- 
though the spatial resolution was lower than the size of individ- 
ual supergranules. They found two regimes taking place in the 
rotation of the photosphere - the quiet Sun and active regions. 
The results showed that the regions occupied by the magnetic 
field display a slower rotation than the non-magnetic vicinity. In 
general, the results of t his study are in agre ement with the re- 
sults of the current one. ISnodgrass & Ulrichl found the mean ro- 
tational rate of magnetic regions with value of (1864+1) ms _1 , 
which, within the statistical errors, agrees with our findings. 
However, the syn odical equatorial rota tion of all Doppler struc- 
tures measured bv lSnodgrass & Ulrichl is (1924+6) m s , which 
is more than 1 % faster than the average rate measured in the 
present study. We could probably explain such a disagreement 
by the different resolution of the used data. In the study of 
ISnodgrass & Ulrichl (1 19901) . the full solar disc was sampled in a 
34x34 pixels array, while in the current one, the size of the solar 
disc was 1000x1000 pixels. We assume that this faster rotational 
rate is due to the effect of undersampling of Doppler structures 
in quiet regions. The results of our study should not be influ- 
enced by this effect, since the supergranular structures are well 
resolved in the data. 

4. Conclusions 

We have verified that the method developed and tested using the 
synthetic data (Paper I) is suitable for application to real data 
obtained by the MDI on-board SoHO and maybe also to the data 
that will be produced by its successor Helioseismic Michelson 
Imager (HMI) on-board the Solar Dynamic Observatory (SDO). 
HMI will have a greater resolution and will cover larger time 
span than two months each year. We verified that the long-term 
evolution of the horizontal velocity fields measured using our 
method is in agreement with generally accepted properties. 

During the periodic analysis of the equatorial area we found 
two suspicious periods in the real data, which are not present in 
the control data set containing the inclination of the solar axis to- 
wards the observer, the quantity that can bias systematically and 
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periodically the results by a few ms _I . The periods of 1.8 year 
and 4.7 years need to be confirmed using a more homogenous 
data set. 

We also found that the presence of the local magnetic field 
generally speeds-up the region occupied by the magnetic field. 
However, we cannot conclude that there exists a dependence 
of this behaviour for different types of sunspots. We can gen- 
erally say that the more evolved types of active regions rotate 
slower than the young ones, however the variance of the typi- 
cal rotation rate is much larger than the differences between the 
rates for each type. We have found that the distribution of ac- 
tive regions rotation is bimodal. The faster-rotating cases corre- 
spond to new and growing active regions. Their almost constant 
rotation speed suggests that they emerge from the base of the 
surface radial shear at 0.95 R e . The decaying and recurrent re- 
gions rotate slower with a wider scatter in their velocities. This 
behaviour suggests that during the sunspot evolution, sunspots 
loose the connection to their magnetic roots. Both regimes al- 
ternate with a period of approximately one Carrington rotation 
in years 2001 and 2002, which suggests that new active regions 
emerge in groups and may have a linked evolution. 
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